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Abstract 

We investigate the relic density n x of non-relativistic long-lived or sta- 
ble particles x i n non-standard cosmological scenarios. We calculate the 
relic abundance starting from arbitrary initial temperatures of the radiation- 
dominated epoch, and derive the lower bound on the initial temperature 
To > m x /23, assuming that thermally produced x particles account for the 
dark matter energy density in the universe; this bound holds for all x annihi- 
lation cross sections. We also investigate cosmological scenarios with modified 
expansion rate. Even in this case an approximate formula similar to the stan- 
dard one is capable of predicting the final relic abundance correctly. Choosing 
the x annihilation cross section such that the observed cold dark matter abun- 
dance is reproduced in standard cosmology, we constrain possible modifications 
of the expansion rate at T ~ m x /20, well before Big Bang Nucleosynthesis. 
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1 Introduction 



One of the most notable recent developments in cosmology is the precise determina- 
tion of cosmological parameters from observations of the large-scale structure of the 
universe, most notably by the Wilkinson Microwave Anisotropy Probe (WMAP). In 
particular, the accurate determination of the non-baryonic cold Dark Matter (DM) 
density [lj, 

0.08 < ttcDMh 2 < 0.12 (95% C.L.) , (1) 

has great influence on particle physics models which possess dark matter candidates 
[21 [3j. The requirement that the predicted DM density falls in the range ([I]) is a 
powerful tool for discriminating between various models and for constraining the 
parameter space of surviving models. 

Many dark matter candidate particles have been proposed. In particular long- 
lived or stable weakly interacting massive particles (WIMPs) with weak-scale masses 
are excellent candidates. In standard cosmology WIMPs decoupled from the thermal 
background during the radiation-dominated epoch after inflation. In this framework 
convenient and accurate analytic approximate solutions for the relic abundance have 
been derived [U [5] . One of the best motivated candidates for WIMPs is the lightest 
neutralino in supersymmetric (SUSY) models. Assuming that the neutralino is the 
lightest supersymmetric particle (LSP) stabilized due to R-parity, its relic abundance 
has been extensively discussed [3J. Similar analyses have also been performed for 
other WIMPs whose existence is postulated in other extensions of the Standard 
Model (SM) of particle physics. In many cases the cosmologically favored parameter 
space of WIMP models can be directly tested at the CERN Large Hadron Collider 
(LHC) in a few years The same parameter space often also leads to rates of 
WIMP interactions with matter within the sensitivity of near-future direct DM 
detection experiments. 

This discussion shows that we are now entering an interesting time where the 
standard cosmological scenario can be examined by experiments at high-energy col- 
liders as well as DM searches [?]. In this respect we should emphasize that the 
relic abundance of thermally produced WIMPs depends not only on their annihila- 
tion cross section, which can be determined by particle physics experiments, but in 
general is also very sensitive to cosmological parameters during the era of WIMP 
production and annihilation. Of particular importance are the initial temperature 
Tq at which WIMPs began to be thermally produced, and the expansion rate of the 
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universe H. 

In the standard cosmological scenario, the expansion rate is uniquely determined 
through the Friedmann equation of general relativity. In this scenario the density of 
WIMPs with mass m x followed its equilibrium value until the freeze-out temperature 
Tp ~ m x /20. Below Tp, interactions of WIMPs are decoupled, and thus the present 
density is independent of T as long as T > T F Q 

It should be noted that in non-standard scenarios the relic density can be larger 
or smaller than the value in the standard scenario. One example is the case where 
To is smaller than or comparable to Tp, which can be realized in inflationary models 
with low reheat temperature. Since in many models the inflationary energy scale 
must be much higher than m x in order to correctly predict the density perturbations 
[8], the standard assumption T > Tp is not unreasonable. On the other hand, the 
constraint on the reheat temperature from Big Bang Nucleosynthesis (BBN) is as low 
as To > MeV [HI [10]. From the purely phenomenological viewpoint, it is therefore 
also interesting to investigate the production of WIMPs in low reheat temperature 
scenarios QU QJl [131 CH] ■ 

The standard scenario also assumes that entropy per comoving volume is con- 
served for all temperatures T < Tp. Late entropy production can dilute the predicted 
relic density [151 [IB]- The reason is that the usual calculation actually predicts the 
ratio of the WIMP number density to the entropy density. On the other hand, if 
late decays of a heavier particle non-thermally produce WIMPs in addition to the 
usual thermal production mechanism, the resulting increase of the WIMP density 
competes with the dilution caused by the decay of this particle into radiation, which 
increases the entropy density [171 [T8l [T9l [201 [TBI [2T] . 

Another example of a non-standard cosmology changing the WIMP relic density 
is a modified expansion rate of the universe. This might be induced by an anisotropic 
expansion [16], by a modification of general relativity [El [22], by additional contri- 
butions to the total energy density from quintessence [23], by branes in a warped 
geometry [21], or by a superstring dilaton [25] . 

These examples show that, once the WIMP annihilation cross section is fixed, 
with the help of precise measurements of the cold dark matter density we can probe 
the very early stage of the universe at temperatures of O(m x /20) ~ 10 GeV. This 
is reminiscent of constraining the early evolution of the universe at T = (9(100) keV 

*Note that Tf can be formally denned in the standard way even if To <Tf- In this case WIMPs 
never were in full equilibrium, and correspondingly never "froze out". 
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using the primordial abundances of the light elements produced by BBN. 

The goal of this paper is to investigate to what extent the constraint (pQ) on 
the WIMP relic abundance might allow us to derive quantitative constraints on 
modifications of standard cosmology. So far the history of the universe has been 
established by cosmological observations as far back as the BBN era. In this paper 
we try to derive bounds on cosmological parameters relevant to the era before BBN. 
Rather than studying specific extensions of the standard cosmological scenario, we 
simply parameterize deviations from the standard scenario, and attempt to derive 
constraints on these new parameters. Since we only have the single constraint (0Q), 
for the most part we only allow a single quantity to differ from its standard value. 
We expect that varying two quantities simultaneously will allow to get the right 
relic density for almost any WIMP annihilation cross section. This has been shown 
explicitly in [T3] for the case that both late entropy production and non-thermal 
WIMP production are considered, even if both originate from the late decay of a 
single scalar field. 

We first analyze the dependence of the WIMP abundance on the initial temper- 
ature Tq of the conventional radiation-dominated epoch. We showed in [14J that for 
fixed T the predicted WIMP relic density reaches a maximum as the annihilation 
cross section is varied from very small to very large values. A small annihilation cross 
section corresponds to a large Tp > T ; in this case the relic density increases with 
the annihilation cross section, since WIMP production from the thermal plasma is 
more important than WIMP annihilation. On the other hand, increasing this cross 
section reduces Tp; once Tp < T a further increase of the cross section leads to 
smaller relic densities since in this case WIMPs continue to annihilate even after the 
temperature is too low for WIMP production. Here we turn this argument around, 
and derive the lower bound on T > m x /23 under the assumption that all WIMPs 
are produced thermally. Note that we do not need to know the WIMP annihilation 
cross section to derive this bound. 

We then examine the dependence of the predicted WIMP relic density on the 
expansion rate in the epoch prior to BBN, where we allow the Hubble parameter to 
depart from the standard value. The standard method of calculating the thermal 
relic density [21 H] is found to be still applicable in this case. Our working hypothesis 
here is that the standard prediction for the Hubble expansion rate is essentially 
correct, i.e. that the true expansion rate differs by at most a factor of a few from 
the standard prediction. We then simply employ a generic Taylor expansion for the 
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temperature dependence of this modification factor; note that the success of standard 
BBN indicates that this factor cannot deviate by more than ~ 20% from unity at 
low temperatures, T < 1 MeV. Similarly, we assume that the WIMP annihilation 
cross section has been determined (from experiments at particle colliders) to have 
the value required in standard cosmology. Our approach is thus quite different from 
that taken in [TJ, where present upper bounds on the fluxes of WIMP annihilation 
products are used to place upper bounds on the Hubble expansion rate during WIMP 
decoupling. The advantage of their approach is that no prior assumption on the 
WIMP annihilation cross section needs to be made, whereas we assume a cross 
section that reproduces the correct relic density in the standard scenario. On the 
other hand, the bounds derived in refs.[7J are still quite weak, allowing the Hubble 
parameter to exceed its standard prediction by a factor > 30; moreover, no lower 
bound on H can be derived in this fashion. 

The remainder of this paper is organized as follows: In Sec. 2 we will briefly review 
the calculation of the WIMP relic abundance assuming a conventional radiation- 
dominated universe, and derive the lower bound on the initial temperature To. In 
Sec. 3 we discuss the case where the pre-BBN expansion rate is allowed to depart 
from the standard one. Using approximate analytic formulae for the predicted WIMP 
relic density for this scenario, we derive constraints on the early expansion parameter. 
Finally, Sec. 4 is devoted to summary and conclusions. 



2 Relic Abundance in the Radiation— Dominated 
Universe 

We start the discussion of the relic density n x of stable or long-lived particles x by 
reviewing the structure of the Boltzmann equation which describes their creation and 
annihilation. The goal of this Section is to find the lowest possible initial temperature 
of the radiation-dominated universe, assuming that the present relic abundance of 
cold dark matter is entirely due to thermally produced \ particles. 

As usual, we will assume that \ is self-conjugatej], \ = x, and that some symme- 
try, for example R-parity, forbids decays of x into SM particles; the same symmetry 
then also forbids single production of x from the thermal background. However, 



^The case x ^ X differs in a non-trivial way only in the presence of a % — \ asymmetry, i.e. if 



4 



the creation and annihilation of x pairs remains allowed. The time evolution of the 
number density n x of particles x i n the expanding universe is then described by the 
Boltzmann equation [2], 

^2L + 3Hn x = -(av)(nl-nlJ, (2) 

where n Xfiq is the equilibrium number density of %, and (av) is the thermally averaged 
annihilation cross section multiplied with the relative velocity of the two annihilating 
X particles. Finally, the Hubble parameter H = R/R is the expansion rate of the 
universe, R being the scale factor in the Friedmann-Robertson- Walker metric. The 
first (second) term on the right-hand side of Eq.([2]) describes the decrease (increase) 
of the number density due to annihilation into (production from) lighter particles. 
Eq.(j2]) assumes that x is m kinetic equilibrium with standard model particles. 

It is useful to rewrite Eq.(T5]) in terms of the scaled inverse temperature x = m x /T 
as well as the dimensionless quantities Y x = n x j s and Y xm = n xeq /s. The entropy 
density is given by s = (27r 2 /45)(y' ;1 , s T 3 , where 

i=bosons i=fermions 

Here gi denotes the number of intrinsic degrees of freedom for particle species i 
(e.g. due to spin and color), and Tj is the temperature of species i. Assuming that 
the universe expands adiabatically, the entropy per comoving volume, sR 3 , remains 
constant, which implies 

f t +3Hs = 0. (4) 
The time dependence of the temperature is then given by 

dx Hx 



dt 



^ x dg* 



(5) 



3#* s dx 

Therefore the Boltzmann equation (T5]) can be written as 

dY x _ (av)s f x dg„ \ 2 _ 2 

dx Hx V 3</„ dx J 1 x x ' eq) ' { ! 

Thermal production of WIMPs takes place during the radiation-dominated epoch, 
when the expansion rate is given by 



with Mpi = 2.4 x 10 GeV being the reduced Planck mass and 

j=bosons i=fermions 

In the following we use H st to denote the standard expansion rate d7j). If the post- 
inflationary reheat temperature was sufficiently high, WIMPs reached full thermal 
equilibrium. This remains true for temperatures well below m x . We can therefore 
use the non-relativistic expression for the \ equilibrium number density, 

n x , eq = 9x (^y 2 e-^ T . (9) 

In the absence of non-thermal production mechanisms, n x < n xm at early times. 
The annihilation rate V = n x (av) then depends exponentially on T, and thus drops 
more rapidly with decreasing temperature than the expansion rate H st of Eq.(j7J) 
does. When the annihilation rate falls below the expansion rate, the number density 
of WIMPs ceases to follow its equilibrium value and is frozen out. 

For T <C m x the annihilation cross section can often (but not always [5]) be 
approximated by a non-relativistic expansion in powers of v 2 . Its thermal average 
is then given by 

(av) = a + b(v 2 ) + 0((v A )) = a + - + ( \) . (10) 

x \x z J 

In this standard scenario, the following approximate formula has been shown [HOE] 
to accurately reproduce the exact (numerically calculated) relic density: 

Y x>00 = Y x (x — >• oo) ~ 1 , (11) 

1.3 m x Mpi X /g :¥ (x F )(a/x F + 3b/x 2 F ) 

with x F = m x /Tp, Tp being the decoupling temperature. For WIMPs, xp ~ 22. 
Here we assume g* ~ g* s and dg*/dx ~ 0. It is useful to express the x mass density 
as fl x = p x j p c , p c = 3HqMp { being the critical density of the universe. The present 
relic mass density is then given by p x = m x n XtOQ = m x soY Xt00 ; here so — 2900 cm -3 
is the present entropy density. Eq. (JTTl) then leads to 

n h 2 - 2 7 x 10 10 Y ( m * ) ~ 8 ' 5 X 1Q ' n XF GeV ' 2 (u) 

where h ~ 0.7 is the scaled Hubble constant in units of 100 km sec" 1 Mpc -1 . We 
defer further discussions of this expression to Sec. 3, where scenarios with modified 
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expansion rate are analyzed. Note that in the standard scenario leading to Eq. ( 112ft . 
the present x re li c density is inversely proportional to its annihilation cross section 
and has no dependence on the reheat temperature. Recall that this result depends 
on the assumption that the highest temperature in the post-inflationary radiation 
dominated epoch, which we denote by To, exceeded Tp significantly. 

On the other hand, if T was too low to fully thermalize WIMPs, the final result 
for Q x will depend on To. In particular, if WIMPs were thermally produced in a com- 
pletely out-of-equilibrium manner starting from vanishing initial abundance during 
the radiation-dominated era, such that WIMP annihilation remains negligible, the 
present relic abundance is given by [1] 



66" 

Y (x -> oo) ~ 0.014 g 2 n; 3 / 2 m x M Pl e- 2xo Xo (« + -)• (13) 



x , 

Note that the final abundance depends exponentially on To, and increases with 
increasing cross section. 

In in-between cases where WIMPs are not completely thermalized but WIMP 
annihilation can no longer be neglected, we have shown [13] that re-summing the 
first correction term S enables us to reproduce the full temperature dependence of 
the density of WIMPs: 

Y ^r^wr Y ^- (14) 



Here 5 < describes the annihilation of WIMPs produced according to Eq. ([T3"l) : 

5{x -> oo) ~ -2.5 x 10~ 4 g^g^mlM^e-^Xo (a + ^ (a + ^ . (15) 

Since 5 is proportional to the third power of the cross section, the re-summed ex- 
pression Y lfJ . is inversely proportional to the cross section for large cross section. In 
ref. [H] we have shown that this feature allows the approximation (fl4|) to be smoothly 
matched to the standard result (fl2"l) . Not surprisingly, as long as we only consider 
thermal x production, decreasing T can only reduce the final x re li c density. 

With the help of these results, we can explore the dependence of the x re li c 
density on To as well as on the annihilation cross section. Some results are shown in 
Fig. [TJ where we take (a) a ^ 0, b = 0, and (b) a = 0, b ^ 0. We choose Y x (xq) = 0, 
m x = 100 GeV, g x = 2 and & = 90. 

The results depicted in this Figure can be understood as follows. For small T , 
i.e. large xq, Eq. flTB"]) is valid, leading to a very strong dependence of Q x h 2 on xo- 
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Figure 1: Contour plots of the present relic abundance Q, x h 2 . Here we take (a) a ^ 0, b = 0, and 
(b) a = 0,6 ^ 0. We choose Y x {x ) = 0, m x = 100 GeV, g x = 2, g.« = 90. The shaded region 
corresponds to the WMAP bound on the cold dark matter abundance, 0.08 < ^cdm^ 2 < O- 1 ^ 
(95% C.L.). 



Recall that in this case the relic density is proportional to the cross section. In this 
regime one can reproduce the relic density (CQ) with quite small annihilation cross 
section, a + 6b/ xq < 10~ 9 GeV -2 , for some narrow range of initial temperature, 
xq < 22.5. Note that this allows much smaller annihilation cross sections than the 
standard result, at the cost of a very strong dependence of the final result on the 
initial temperature T . 

In this Section we set out to derive a lower bound on T . In this regard the 
region of parameter space described by Eq. (fl3|) is not optimal. Increasing the x 
annihilation cross section at first allows to obtain the correct relic density for larger 
Xq, i.e. smaller T . However, the correction 5 then quickly increases in size; as noted 
earlier, once \5\ > Yq a further increase of the cross section will lead to a decrease 
of the final relic density. The lower bound on To is therefore saturated if Q x h 2 as a 
function of the cross section reaches a maximum. From FigffJ we read off 

T > m x /23 , (16) 

if we require Q x h 2 to fall in the range ([T]). 

We just saw that in the regime where this bound is saturated, the final relic den- 
sity is (to first order) independent of the annihilation cross section, d(Q x h 2 )/d(av) = 
0. If T is slightly above the absolute lower bound (jTBT) . the correct relic density can 
therefore be obtained for a rather wide range of cross sections. For example, if 



S 



x = 22.5, the entire range 3 x 1(T 10 GeV~ 2 < a < 2 x KT 9 GeV" 2 is allowed. Of 
course, the correct relic density can also be obtained in the standard scenario of 
(arbitrarily) high T , if a + 36/22 falls within ~ 20% of 2 x KT 9 GeV~ 2 . 



3 Relic Abundance for Modified Expansion Rate 

In this section we discuss the calculation of the WIMP relic density n x in modi- 
fied cosmological scenarios where the expansion parameter of the pre-BBN universe 
differed from the standard value H st of Eq.([7]). For the most part we will assume 
that WIMPs have been in full thermal equilibrium. Various cosmological models 
predict a non-standard early expansion history [221 E3J EH EH] . Here we analyze to 
what extent the relic density of WIMP Dark Matter can be used to constrain the 
Hubble parameter during the epoch of WIMP decoupling. As long as we assume 
large T we can use a modification of the standard treatment [U E] to estimate the 
relic density for given annihilation cross section and expansion rate. We will show 
that the resulting approximate solutions again accurately reproduce the numerically 
evaluated relic abundance. 

Let us introduce the modification parameter A(x), which parameterizes the ratio 
of the standard value H st (x) to the assumed H(x): 

(17) 

Note that A > 1 means that the expansion rate is smaller than in standard cos- 
mology. Allowing for this modified expansion rate, the Boltzmann equation ([6]) is 
altered to 

dY x _ 4tt (av)A{x) , 2 2 , 

where 



fgl \ 3g* s dx 

Following refs. [H E] , we can obtain an approximate solution of this equation by 
considering the differential equation for A = Y x — Y xm . For temperatures higher 

be ignored: 



= (l - ^% ) • (19) 



than the decoupling temperature, Y x tracks Y xm very closely and the A 2 -term can 



dA dF cq 4tt w G(x)(av)A(x) . a \ 

z ~ - —=m x M vx K A / ^ ; (2K eq A . 20 

dx dx v/90 x x 2 v x,q ' v ; 



Here dY Xfiq /dx ~ — V XjCq for i> 1. In order to keep |A| small, the derivative dA/dx 
must also be small, which implies 

(21) 



A 



x 



^/V90)m x M Pl G(x)(av}A(x) 
This solution is used down to the freeze-out temperature Tp, defined via 

A{x F ) = £Y x>eq {x F ) , (22) 

where £ is a constant of order of unity. This leads to the following expression: 

^45 



Xp = In 




(av)A(x) ( x dg* s 

£m x M Pl g x 1 - — 

vr 5 Vxgjx) V 3#* s dx 



(23) 



X=Xp 



which can e.g. be solved iteratively. In our numerical calculations we will choose 

f =V2-i mi. 

On the other hand, for low temperatures (T < Tp), the production term oc Y xeq 
in Eq. fll8p can be ignored. In this limit, Y x ~ A, and the solution of Eq. (|T8|) is given 
by 

1 1 4vr 



A(x F ) A(x oo) 
where the annihilation integral is defined as 

G{x)(av)A{x) 



m x Mp\I(x F ) 



I(x F ) 



dx 



Xp 



Assuming A(x — > oo) A(x F ), the final relic abundance is 

1 



Y x ,oo = Y x (x -> oo) 



(24) 



(25) 



(26) 



(4vr/v / 90)m x M P1 J(x F ) 

Plugging in numerical values for the Planck mass and for today's entropy density, 
the present relic density can thus be written as 

8.5 x 10~ u 



Q x h 2 



(27) 



I(xp) GeV 2 ' 

The constraint ([1]) therefore corresponds to the allowed range for the annihilation 
integral 



7.1 x 1(T 1U GeV < I(x F ) < 1.1 x 10" 9 GeV" 



(28) 



The standard formula ( fT2l) for the final relic density is recovered if A(x) is set to 
unity and G(x) is replaced by the constant \/g*(x F ). 
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The further discussion is simplified if we use the normalized temperature z = 
T/m x = 1/x, rather than x. Phenomenologically A(z) can be any function subject 
to the condition that A(z) approaches unity at late times in order not to contradict 
the successful predictions of BBN. We need to know A(z) only for the interval from 
around the freeze-out to BBN: zbbn ~ 10~ 5 — 1CT 4 < z < zp ~ 1/20. This suggests 
a parameterization of A(z) in terms of a power series in (z — z FjSt ): 

A{z) = A(z FjSt ) + (z- z F)St )A'{z F)St ) + -(z - z F , st ) 2 A"{z F , Bt ) , (29) 

where z F>st is the normalized freeze-out temperature in the standard scenario and 
a prime denotes a derivative with respect to z. The ansatz (12"9"|) should be of quite 
general validity, so long as the modification of the expansion rate is relatively modest. 
This suits our purpose, since we wish to find out what constraints can be derived on 
the expansion history if standard cosmology leads to the correct WIMP relic density. 
We further introduce the variable 

k = A(z -> 0) = A(z F)St ) - z F ^A'(z F)St ) + ^z 2 FtSt A"{z F , st ) , (30) 

which describes the modification parameter at late times. Since zbbn is almost 
zero, we treat k as the modification parameter at the era of BBN in this paperlfl 
Deviations from k = 1 are conveniently discussed in terms of the equivalent number 
of light neutrino degrees of freedom N v . BBN permits that the number of neutrinos 
differs from the standard model value N v = 3 by 8N V = 1.5 or so We therefore 
take the uncertainty of k to be 20%. In the following we treat A(z F)St ), A'(z FyS t) and 
k as free parameters; A"(z F ^ st ) is then a derived quantity. 

Note that we allow A(z) to cross unity, i.e. to switch from an expansion that 
is faster than in standard cosmology to a slower expansion or vice versa. This is 
illustrated in Fig. [2], which shows examples of possible evolutions of A(z) as function 
of z for z F = 0.05. Here we take k = 1.2 (left frame) and k = 0.8 (right). In each case 
we consider scenarios with A(z F ) = 1.3 (slower expansion at T F than in standard 
cosmology) as well as A(z F ) = 0.7 (faster expansion); moreover, we allow the change 
of A at z = z F to be either positive or negative. However, we insist that H remains 
positive at all times, i.e. A(z) must not cross zero. This excludes scenarios with 
very large positive A'(z FtSt ), which would lead to A < at some z < z F . Similarly, 

* Presumably the Hubble expansion rate has to approach the standard rate even more closely 
for T < Tbbn- However, since all WIMP annihilation effectively ceased well before the onset of 
BBN, this epoch plays no role in our analysis. 
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0.01 0.02 0.03 0.04 0.05 0.01 0.02 0.03 0.04 0.05 

z = T/m^ z = T/m x 

(a) (b) 

Figure 2: Examples of possible evolutions of the modification parameter A(z) as function of z 
for z F — 0.05. Here we take k = 1.2 (left frame) and k = 0.8 (right). In each frame we choose 
A(z F ) = 1.3, A'(z F ) = -3 (thick line), A(z F ) = 1.3, A'(z F ) = 9 (dashed), = 0.7, A'(z F ) = 

-3 (dotted), A(z F ) = 0.7, A'(z F ) = 9 (dot-dashed). 



demanding that our ansatz (1291) remains valid for some range of temperatures above 
Tp excludes scenarios with very large negative A'(zp }B t). We will come back to this 
point shortly. 

Eq. fTSBl shows that Zp ^ Zp ;St [xp ^ xp >st ) if A(zp) ^ 1. This is illustrated by 
Fig. [31 which shows the difference between xp and xp }S t in the (A(zF )St ), ^'(^st)) 
plane. Here we take parameters such that Q x h 2 = 0.099 in the standard cosmology, 
which is recovered for A{zp iS t) = 1, A'(zp^t) = 0. Due to the logarithmic dependence 
on A, xp (or zp) differs by at most a few percent from its standard value if A{zf^\) is 
(9(1). Since Tp only depends on the expansion rate at Tp, it is essentially insensitive 
to the derivative A'(zp )St ). 

In our treatment the modification of the expansion parameter affects the WIMP 
relic density mostly via the annihilation integral (T25]) . In terms of the normalized 
temperature z, the latter can be rewritten as 

I( ZF ) = j F Az G(z)(av)A(z) . (31) 
Jo 

One advantage of the expansion (f29l) is that this integral can be evaluated analyti- 
cally: 



I(z F ) ~ G{zp) [k{az F + 364) + W^,*) - z F , st A" (z FjBt )) (~z 2 F + 2bz 



A"(z F>at ) fa ,3b 4 

—z% H z% 

3 F 2 F 
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(32) 



< 
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: 2.0*10 9 GeV 2 
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Figure 3: Contour plot of xp — xp^t in the (A(zF,st), A'(zF, s t)) plane. Here we take a = 2.0 x 10~ 9 
GeV -2 , b = 0, m x = 100 GeV, g x = 2, g* = 90 (constant) and k = 1. This parameter set produces 
XF,st = 22.0 and il x h 2 = 0.099 for the standard approximation. 



n z h 2 (approx) / n z h 2 (exact) n z h 2 (approx) / n x h 2 (exact) 




Figure 4: Ratio of the analytic result of the relic density to the exact value in the (A(zF.st), 
A'{zF,st)) plane for a = 2.0 x 10~ 9 GcV~ 2 , 6 = (left frame) and for a = 0, b = 1.5 x 10" 8 GeV" 2 
(right). The other parameters are as in Fig. [3] 



Here we have assumed that G(z) varies only slowly. 

Before proceeding, we first have to convince ourselves that the analytic treatment 
developed in this Section still works for A ^ 1. This is demonstrated by Fig. 0], 
which shows the ratio of the analytic solution obtained from Eqs. (T27|) and (|32|) 
to the exact one, obtained by numerically integrating the Boltzmann equation (|T8l) . 
assuming constant g*. We see that our analytical treatment is accurate to better 
than 1%, and can thus safely be employed in the subsequent analysis. 
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Figure 5: Contour plots of the relic abundance in the (A(zF,nt)i A'{zF,st)) plane. Here we choose 
(a) a = 2.0 x 1CT 9 GeV~ 2 , b = 0, k = 1; (b) a = 0, b = 1.5 x 1CT 8 GeV~ 2 , k = 1; (c) a = 2.0 x 10~ 9 
GeV^ 2 , b = 0, k = 1.2; (d) a = 2.0 x 10~ 9 GeV~ 2 , b = 0, k = 0.8. The other parameters are as in 

Fig. m 



We are now ready to analyze the impact of the modified expansion rate on the 
WIMP relic density. In Fig. [5], we show contour plots of Q x h 2 in the (A(zF, s t), 
A'(zF,st)) plane. Recall that large (small) values of A correspond to a small (large) 
expansion rate. Since a smaller expansion rate allows the WIMPs more time to 
annihilate, A > 1 leads to a reduced WIMP relic density, whereas A < 1 means 
larger relic density, if the cross section is kept fixed. 

However, unlike the freeze-out temperature, the annihilation integral is sensitive 
to A(z) for all z < zp. Note that A'(zF, B t) > implies A(z) < A(zF, s t) for z < 
ZF t st — z f- A positive first derivative, A'(zF >s t) > 0, can therefore to some extent 
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compensate for A(z Fs t) > 1; analogously, a negative first derivative can compensate 
for A(z Fst ) < 1. This explains the slopes of the curves in Fig. [51 Recall also that 
A'(zF tSt ) = does not imply a constant modification factor A(z)\ rather, the term 
oc A"(zF,$t) in Eq. (}2"9]) makes sure that A approaches k as z — > 0. This explains why 
a change of A by some given percentage leads to a smaller relative change of Q x h 2 , as 
can be seen in the Figure. This also illustrates the importance of ensuring appropriate 
(near-standard) expansion rate in the BBN era. Finally, since the expansion rate 
at late times is given by H st /k, bigger (smaller) values of k imply that the x rei i c 
density is reduced (enhanced). 

Fig. [5] shows that we need additional physical constraints if we want to derive 
bounds on A(zF, s t) and/or A'(z F)S t). Once the annihilation cross section is known, 
the requirement (PQ) will single out a region in the space spanned by our three new 
parameters (including k) which describe the non-standard evolution of the universe, 
but this region is not bounded. Such additional constraints can be derived from 
the requirement that the Hubble parameter should remain positive throughout the 
epoch we are considering. As noted earlier, requiring H > for all T < Tp, s t leads 
to an upper bound on A'(z Fst ); explicitly, 

2 (a(z f &) + y/kA(z FjBt )) 
A'(z F ,st) < -± L ■ (33) 

z F) st 

On the other hand, a lower bound on A'(z Fst ) is obtained from the condition that 
the modified Hubble parameter is positive between the highest temperature Tj where 
the ansatz fl29|) holds and T Fst : 

(34) 



A'(z Fst ) > - [ (2 - A(z Fst ) + k (— - - 

_ Zi — Z Fs t \ Z Ft st J \ Z F ,st z i 

for (1 - z F , s t/zi) 2 k < A(z Fet ), and 

2 (A(z F , st ) - y/kA(z FiSb j) 
A'(z F ,st) > L , (35) 

Z F ,St 



for A(z F>st ) < (1 - z F , s t/zi) 2 k, where Zi = Ti/m 

Evidently the lower bound on A^^st) depends on Zi, i.e. on the maximal tem- 
perature where we assume our ansatz (1291) to be valid. In ref. [H] we have shown that 
in standard cosmology (A = 1) essentially full thermalization is already achieved for 
x i ~ x f ~ 0.5, even if n x (xi) = 0. However, it seems reasonable to demand that H 
should remain positive at least up to X{ = x F —(a few). In Fig. [6] we therefore show 
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Figure 6: Contour plots of the relic abundance fl x h 2 in the (A(zF, s t),A'(zF, a t)) plane. The 
dashed line corresponds to the upper bound on A'(zF.st)- The dotted lines correspond to the lower 
bounds calculated for xp, s t — %i = 4, 10. We take a = 2.0 x 10 -9 GeV~ 2 , 6 = (left frame) and 
a = 0, b = 1.5 x 10 -8 GeV~ 2 (right frame). The other parameters are as in Fig [31 



the physical constraints on the modification parameter A(z) for Xp >s t — Xi = 4, 10 
and k — 1. The dashed and dotted lines correspond to the upper and lower bounds 
on A'(z Fst ), described by Eq. (l33l and Eqs. tjM]) . (|35|) . respectively We see that when 
%F,st - ^ = 4 the allowed region is 0.4 < A(z F ,st) ^ 6.5 with -60 < A'(z F>st ) < 400 
for b = (left frame), and 0.4 < A(z F , st ) < 4.5 with -60 < A'(z FjSt ) < 300 for a = 
(right frame). When x Fst — Xi = 10, the lower bounds are altered to 0.6 < A(z Fst ), 
-10 < A'(z F>st ) for b = (left frame), and 0.6 < A{z Fst ), -20 < A'(z Fst ) for a = 
(right frame). Note that the lower bounds on A(z Fst ), which depend only weakly on 
Xi so long as it is not very close to x F , are almost the same in both cases, which also 
lead to very similar relic densities in standard cosmology. However, the two upper 
bounds differ significantly The reason is that large values of A(z Fst ), i.e. a strongly 
suppressed Hubble expansion, require some degree of finetuning: One also has to 
take large positive A'(z Fst ), so that A becomes smaller than one for some range of 
z values below z F , leading to an annihilation integral of similar size as in standard 
cosmology. Since the b— terms show different z F dependence in the annihilation in- 
tegral (I3~2"l) . the required tuning between A(z F>st ) and A'(z F>st ) is somewhat different 
than for the a— terms, leading to a steeper slope of the allowed region. This allowed 
region therefore saturates the upper bound f)33p on the slope for somewhat smaller 
A{z FjBt ). 

The effect of this tuning can be seen by analyzing the special case where A"(z Fst ) = 
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Figure 7: Contour plots of the relic abundance Vl x h 2 in the (A(zp, s t), k) plane for A"(zF, s t) = 0. 
The dotted lines correspond to the lower bounds of A(zp^ t ), calculated for XF,st — x*i — 4, 10. We 
take a = 2.0 x 10~ 9 GeV" 2 , 6 = (left frame) and a = 0, b = 1.5 x 10~ 8 GeV~ 2 (right frame). 
The other parameters are as in Fig[3J 



0. The modification parameter then reads 

A(z) = A{ZF ' st) ~ k z + k . (36) 

Z F,Bt 

Note that A is now a monotonic function of z, making large cancellations in the 
annihilation integral impossible. Imposing that A(z) remains positive for zp, s t < 
z < Zi leads to the lower limit 

A(z F>st ) > (l - ^) k . (37) 

There is no upper bound, since A(z) is now automatically positive for all z G [0, zp^ st } 
if A(zF,st) an d ^4(0) = k are both positive. Fig. [7] shows constraints on the relic 
abundance in the (A(zp )St ), k) plane for A"(zF,st) — 0- The dotted lines correspond 
to the lower bounds (1371) on A(zF,st) f° r x F,st — =4, 10. As noted earlier, k 
is constrained by the BBN bound. This leads to the bounds 0.5 < A(z F , s t) i$ 1-8 
for 6 = (left frame), and 0.65 < A(zF, s t) ~ 1-6 for a = (right frame), when 
%F,st — Xi = 10. Evidently the constraints now only depend weakly on whether the 
a— or b— term dominates in the annihilation cross section. As the initial temperature 
is lowered, the impact of the constraint (|37|) disappears. 

So far we have assumed in this Section that the reheat temperature was high 
enough for WIMPs to have attained full thermal equilibrium. If this was not the 
case, the initial temperature as well as the suppression parameter affects the final 
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Figure 8: Contour plot of the relic abundance Q, x h 2 in the (A(zF,st), xo) plane. Here we choose 
a = 2.0 x 1CT 9 GeV~ 2 , b = 0, k = l,A"(z FjSt ) = 0. The other parameters are as in Fig. [3 
The shaded region corresponds to the WMAP bound on the cold dark matter abundance, 0.08 < 
Ocdm/i 2 < 0.12 (95% C.L.). 



relic abundance. Here we show that the lower bound on the reheat temperature 
derived in the previous Section survives even in scenarios with altered expansion 
history as long as WIMPs were only produced thermally 

This can be understood from the observation that the Boltzmann equation with 
modified expansion rate is obtained by replacing (av ) in the radiation-dominated 
case by (av)A. Increasing (decreasing) A therefore has the same effect as an in- 
crease (decrease) of the annihilation cross section. Since the lower bound on To was 
independent of a (more exactly: we quoted the absolute minimum, for the optimal 
choice of a), we expect it to survive even if A(z) ^ 1 is introduced. 

This is borne out by Fig. [HI which shows the relic abundance Q x h 2 in the (A(zF,st), 
xq) plane for the simplified case A"{zF^t) = 0; similar results can be obtained for 
the more general ansatz (1291) . The shaded region corresponds to the bound ([I]) on 
the cold dark matter abundance. As expected, this figure looks similar to Fig. [T] if 
the annihilation cross section in Fig. [1] is replaced by ^(^st). The maximal value of 
xq consistent with the WMAP data remains around 23 even in these scenarios with 
modified expansion rate. Fig. [HJ also shows that A(zF, s t) 1 is allowed for some 
narrow range of initial temperature T < Tp. This is analogous to the low cross 
section branch in Fig. [TJ 
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4 Summary and Conclusions 



In this paper we have investigated the relic abundance of WIMPs x, which are non- 
relativistic long-lived or stable particles, in non-standard cosmological scenarios. 
One motivation for studying such scenarios is that they allow to reproduce the ob- 
served Dark Matter density for a large range of WIMP annihilation cross sections. 
Our motivation was the opposite: we wanted to quantify the constraints that can 
be obtained on parameters describing the early universe, under the assumption that 
thermally produced WIMPs form all Dark Matter. Wherever necessary, we fixed 
particle physics quantities such that standard cosmology yields the correct relic den- 
sity. 

Specifically, we first considered scenarios with low post-inflationary reheat tem- 
perature, while keeping all other features of standard cosmology (known particle 
content and Hubble expansion parameter during WIMP decoupling; no late entropy 
production; no non-thermal WIMP production channels). If the temperature was so 
low that WIMPs could not achieve full thermal equilibrium, the dependence of the 
abundance on the mass and annihilation cross section of the WIMPs is completely 
different from that in the standard thermal WIMP scenario. In particular, if the 
maximal temperature To is much less than the decoupling temperature Tp, n x re- 
mains exponentially suppressed. By applying the observed cosmological amount of 
cold dark matter to the predicted WIMP abundance, we therefore found the lower 
bound of the initial temperature To > m x /23. One might naively think that this 
bound could be evaded by choosing a sufficiently large WIMP production (or anni- 
hilation) cross section. However, increasing this cross section also reduces T F . For 
sufficiently large cross section one therefore has T F < T again; in this regime the 
relic density drops with increasing cross section. Our lower bound is the minimal 
T required for any cross section; once the latter is known, the bound on T might 
be slightly sharpened. As a by-product, we also noted that the final relic density 
depends only weakly on the annihilation cross section if T is slightly above this 
lower bound. 

We also investigated the effect of a non-standard expansion rate of the universe 
on the WIMP relic abundance. In general the abundance of thermal relics depends 
on the ratio of the annihilation cross section to the expansion rate; the latter is 
determined unambiguously in standard cosmology. We found that even for non- 
standard Hubble parameter the relic abundance can be calculated accurately in 
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terms of an annihilation integral, very similar to the case of standard cosmology. We 
assumed that the WIMP annihilation cross section is such that the standard scenario 
yields the observed relic density, and parameterized the modification of the Hubble 
parameter as a quadratic function of the temperature. In this analysis it is crucial to 
make sure that at low temperatures the Hubble parameter approaches its standard 
value to within ~ 20%, as required for the success of Big Bang Nucleosynthesis 
(BBN). 

Keeping the annihilation cross section fixed and allowing a 20% variation in the 
relic density, roughly corresponding to the present "2<t" band, we found that the 
expansion of the universe at T = Tp might have been more than two times faster, or 
more than six times slower, than in standard cosmology. These large variations of 
H(Tp) can only be realized by finetuning of the parameters describing H(T < Tp). 
However, even if we forbid such finetuning by choosing a linear parameterization for 
the modification of the expansion rate, a 20% variation of Q x h 2 allows a difference 
between H(Tp) and its standard expectation of more than 50%. This relatively weak 
sensitivity of the predicted fl x h 2 on H{Tp) is due to the fact that the relic density 
depends on all H(T < T F ); as stressed above, we have to require that H(T <C T F ) 
approaches its standard value to within ~ 20%. The fact that determining Q x h 2 will 
yield relatively poor bounds on H(Tp) remains true even if the annihilation cross 
section is such that a non-standard behavior of H(T) is required for obtaining the 
correct x renc density. Finally, we showed that the absolute lower bound on the 
temperature required for thermal x production is unaltered by allowing H(T) to 
differ from its standard value. 

Of course, in order to draw the conclusions derived in this article, we need to 
convince ourselves that WIMPs do indeed form (nearly) all Dark Matter. This 
requires not only the detection of WIMPs, e.g. in direct search experiments; we 
also need to show that their density is in accord with the local Dark Matter density 
derived from astronomical observations. To that end, the cross sections appearing in 
the calculation of the detection rate need to be known independently. This can only 
be done in the framework of a definite theory, using data from collider experiments. 
For example, in order to determine the cross section for the direct detection of 
supersymmetric WIMPs, one needs to know the parameters of the supersymmetric 
neutralino, Higgs and squark sectors [3]. We also saw that inferences about H(T F ) 
can only be made if the WIMP annihilation cross section is known. This again 
requires highly non-trivial analyses of collider data [27], as well as a consistent 
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overall theory. We thus see that the interplay of accurate cosmological data with 
results obtained from dark matter detections and collider experiments can give us 
insight into the pre-BBN universe, which to date remains unexplored territory. 
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